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ABSTRACT 

SN1957D, located in one of the spiral arms of M83, is one of the small number 
of extragalactic supernovae that has remained detectable at radio and optical 
wavelengths during the decades after its explosion. Here we report the first 
detection of SN1957D in X-rays, as part of a 729 ks observation of M83 with 
Chandra. The X-ray luminosity (0.3 - 8 keV) is 1-71 2 ,; 4 x 10 37 ergss -1 . The 
spectrum is hard and highly self-absorbed compared to most sources in M83 and 
to other young supernova remnants, suggesting that the system is dominated at 
X-ray wavelengths by an energetic pulsar and its pulsar wind nebula. The high 
column density may be due to absorption within the SN ejecta. HST WFC3 
images resolve the supernova remnant from the surrounding emission and the 
local star field. Photometry of stars around SN1957D, using WFC3 images, 
indicates an age of less than 10 7 years and a main sequence turnoff mass more 
than 17 M . New spectra obtained with Gemini-South show that the optical 
spectrum continues to be dominated by broad [O III] emission lines, the signature 
of fast-moving SN ejecta. The width of the broad lines has remained about 
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2700 kms" 1 (FWHM). The [O III] flux dropped precipitously between 1989 
and 1991, but continued monitoring shows the flux has been almost constant 
since. In contrast, radio observations over the period 1990-2011 show a decline 
rate S u ~ £~ 4 - , far steeper than the rate observed earlier, suggesting that the 
primary shock has overrun the edge of a pre-SN wind. 

Subject Headings: galaxies: individual (M83) - galaxies: ISM - radio contin- 
uum: galaxies - supernova remnants - X-rays: general - X-rays: individual (SN 
1957D, M83) - supernovae: individual (SN 1957D) 



1. Introduction 



M83 (NGC 5236) is a nearby grand-design spiral with prolific star formation that 
has been the host to six historical supernovae (SNe), including SN1957D. Relatively lit- 
tle is known about the SN1957D. It was observed at a photographic magnitude m pg < 15 



( iKowal fc Sargentlll97ll ). which corresp onds to an absolu te magnitude of -13.3 using an as- 
sumed distance to M83 of 4.61 Mpc ( ISaha et al.l 120061 ). It was probably well past peak 
brightness at the time when it was observed. The SN is located at the inner edge of a 
prominent northern spiral arm about 3 kpc from the nucleus of the galaxy. No spectra 
during outburst are known, and so we do not know the SN type from original observations. 
However, because of its proximity to active star formation, a core-collapse SN of some type 
is strongly suspected. 

SN1957D was first recovered as a young supernova remnant (SN R) at radio wavelengths 
in 1981, 24 years after the SN explosion, by lCowan fc Branch! ( 119851 ) . A nonthermal source, 
it had a monochromatic radio luminosity similar to that of Cas A at that time. Its flux 
declined over t he next few years, as do most old SNe/young SNRs that are detected in 
this age range Jstockdale et al.ll2006h . Of the six historical SNe in M83, three - SN1923A, 
1950B, and 1957D - remained detectable as radio sources years after the SN event. SN1983N, 



1 Based on observations made with NASA's Chandra X-ray Observatory, the NASA/ESA Hubble Space 
Telescope, the 6.5 meter Magellan Telescopes located at Las Campanas Observatory, and the Gemini Obser- 
vatory. NASA's Chandra Observatory is operated by Smithsonian Astrophysical Observatory under contract 
# NAS83060 and the data were obtained through program GOl-12115. The HST observations were obtained 
at the Space Telescope Science Institute, which is operated by the Association of Universities for Research 
in Astronomy, Inc. (AURA), under NASA contract NAS 5-26555. The new HST observations were obtained 
through programs # GO-12513 and GO-12683. Data in the HST archive from program GO-11360 was also 
used. The ground-based observations were obtained through NOAO which is operated by Association of 
Universities for Research in Astronomy, Inc. for the National Science Foundation. 
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classified as a type lb SN, was detect ed at radio wavelengths just after the expl osion, but 
faded below detectability very quickly (jCowan Sz Branchlll985l ; lMaddox et al. hood ). The two 
other historical SNe, SN1945B and SN1968L, have not been detected in the radio. SN1968L 
lies is in a region of bright, diffuse radio emission near the nucleus, making detection of a 
radio SNR much more difficult than the others with arcsecond scale angular resolution^ 

SN1957D was rec overed at optical wavelengths in 1987, 30 years after the SN explosion, 
by lLong et al.l (119891 ) as an [O III] emis sion- line source. By 1991, its flux declined by a 
factor of five when it was reobserved by lLong et al.l (119921 ). Spectra obtained during this 
period showed emission in [O II] AA3727, 3729, [O III] AA4959, 5007 and [O I] AA6300, 6363 
but little else. Line widths were 2400-2900 km s~ l (FWHM) when the signal-to-noise ratio 



was h igh enough to measure the line width (lLong et al.lll989l ; iTuratto et al.lll989l ; lLong et al. 



19921 ). The spectra resemble those of other young SNRs thought to have arisen from core- 
collapse SNe, such as Cas A and G292+1.8 in the Galaxy and N132D, and E0102-72.3 
in the Magellanic Clouds. Emission in these O-rich SNRs is usually interpreted as aris- 
i ng from shock-heated knot s of metal-enriched ejecta from the explosion of a massive star 
JSutherland fc Dopital[l99~5T ). 



SN1957D h as not been detecte d previously at X-ray wavelengths, and in particular it was 
not detected by ISoria &: Wul (120031 ) who analyzed a 50 ks observation of M83 obtained with 
Chandra in 2000. However, it is clearly detected as a very hard X-ray point source in new, 
much deeper, Chandra observations with a total exposure time of 729 ks. This detection, 
along with new radio and optical observations of the SNR and its local environment, are the 
subject of this report. In the next three sections, we describe the results of our analysis of 
the observations at X-ray, radio, and optical wavelengths. Then in Section [5], we compare 
SN1957D to other very young SNRs in order to obtain an overall picture of the SNR. Finally, 
in Section [6J we summarize the main results. 



2. X-ray Observations - Reductions and Analysis 

The new Chandra observations of M83 were obtained in ten sections, starting in late 
2010 and continuing throughout 2011 as part of a large program (12620596; Long P.I.) to 
characterize SNRs, point sources and diffuse emission in M83 to a limiting luminosity of 
about 5 x 10 35 ergs s _1 (for spectral parameters typical of X-ray bi naries). The com bined 



exposure totals 729 ks. A detailed observing log appears in Table 1 of ISoria et al.l (120121 ) . We 



used the back-illuminated S3 chip in order to maximize the response at soft X-ray energies 



2 SN1968L has been tentatively identified optically in HST/WFC3 images jPopita et al.lboiol) . 
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and we obtained the data in the "very fa int" mode to optimize background subtraction. We 



have used tools within CIAO Version 4.4 (IFruscione et al.ll2006l ) to reduce the data, especially 



to create response and a rea corrected ancillary files needed for spectral fitting and we have 



used XSPEC Version 12 (jArnaudlll996l ) to characterize the X-ray spectrum 



A 3-color X-ray image of the region of M83 containing SN1957D is shown in the upper 
left panel of Fig. [TJ For comparison, an optical emission line image obtained at the 6.5m 
Magellan telescope at Las Campanas Observatory and optical continuum and emission line 
images obtained with HST are also shown. There is a hard (blue) X-ray source coincident 
with an [O III]-dominated, point-like emission nebula that we identify as the young SNR 
from SN1957D (see following sections). The SNR is also detected as a non-thermal radio 
source with ATCA (see Section [3]). 

The uncertainty in the absolute astrometry of pipeline-processed Chandra images is 0'.'6 
(90% confidence) within 3' of the aim point Jf] The statistical uncertainty in the position of 
the source simply due to the size of the point spread function and number of X-ray events 
detected is of order Qf!l (la). Our initial position estimate for the X-ray source identified at 
SN1957D was about Qf!3 from the ATCA radio position for SN1957D, so within the errors 
the positions are identical. However, the absolute astrometry of the Chandra images can be 
improved by using the ATCA radio data. Of the 65 sources observed with ATCA at 5 GHz, 
25 lie within 1" of sources in our preliminary version of the M83 X-ray point source catalog. 
If we exclude SN1957D and correct the astrometry of the Chandra positions by the average 
position offset of the 24 remaining sources, the difference between the corrected X-ray source 
and radio position of SN1957D is 0'.'14, is very close to the statistical error. The corrected 
(J2000) position for the object we identify with SN1957D is: 

X-ray : 13 h 37 m 03.584 s ; -29 o 49'40.9i" 

To obtain the X-ray spectrum of SN1957D, we defined a circular source region with 
radius 2'.'5, and a background annulus between radii 3" and 6'.'5. We extracted and combined 
source and background spectra from each observation, together with their response and an- 
cillary response matrices. We have fit the spectrum, which is shown in Fig. [2] and contains 
140 counts after background subtraction, to power-law and thermal models. Due to the 
limited number of counts, we have used c-statistics rather tha n y 2 statistics in the m odel fit- 



ting. Fixing foreground Galactic absorption at 4 x 10 cm ( IKalberla et al.ll2005l ). we find 



J See the analysis provided by the CXC at http://cxc.harvard.edu/cal/ASPECT/celmon/ 



-5 - 



a power law with a photon power law index 7 of 1.41J g and iV#(M83) of 2.0l^ x 10 22 cm- 2 @ 
The derived 0.3-8 keV luminosity is I.7I03 x er S s s_1 - The spectrum can also be mod- 
eled equally well as with a thermal model, but in this case we obtain only a lower limit of 
10 keV to the temperature. If X-ray emission from SN1957D is from a thermal plasma with 
this effective temperature, we would not expect to see lines in spectra with so few counts. 

We searched for evidence of a source at the position of M83 in archival Chandra data 
obtained in 2000 and 2001 with a total exposure of 60,800 s. The source was not detected 
although there was a slight excess at the known position of the SN. We estimate the 0.3-8 keV 
count rate to have been 9.4 ± 7.9 x 10~ 5 counts s _1 , corresponding to a 3 a upper limit of 
3.3 x 10~ 4 counts s -1 . In 2010, the count rate was 1.8 ± 0.2 x 10~ 4 counts s" 1 . As a result, 
we cannot exlude the possibility that the source was as bright in the archival data as it is 
today, but we can rule out the possibility that it was much brighter in 2000-2001. 



3. Radio Observations — Reductions and Analysis 



M83 was observed with the Australia Telescope Compact Array (ATCA) in the 6A 
configuration on 2011 April 28-30. On each of the three days, a full 12-hour synthesis was 
obtained, simultaneously recording 2 GHz bandwidths centered on 5.5 GHz and 9 GHz in 
all four polarization products. Regular scans on the nearby phase calibrator PKS 1313- 
333 were scheduled throughout the observations, as well as scans on the ATCA primary 
flux calibrator PKS 1934-638. Standard c alibration and edi ting procedures were carried out 
using the MIRIAD data analysis package (ISault et al.lll995l ). Following calibration, the data 
were exported for imaging, self-calibration and prima ry beam correction usi ng the Common 
Astronomy Software Applications (CASA) package (iMcMullin et al.l 120071 ). The 5.5 GHz 
data were imaged with uniform weighting and had angular resolution 2.52"x 1.06" at a 
position angle of 179°. The 9 GHz map was imaged with Briggs weighting (Robust = 0) 
and convolved to match the 5.5 GHz resolution. Coordinates were obtained by fitting an 
elliptical Gaussian to the pixels with values greater than half the peak, using the CASA task 
IMFIT. We detected about 64 and 23 sources at 5.5 and 9 GHz respectively, including, at 
both frequencies, SN1957D. A complete discussion of the ATCA radio observations of M83 
will appear elsewhere. 



SN1957D is unresolved in both the 5.5 GHz and 9 GHz images and has a measured 



4 All of the errors quoted here represent 90% confidence intervals. 
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position (J2000) of 

Radio : 13 h 37 m 03.586 s ± 0.003 s ; -29°49'41.0" ± 0.1'.' 

Flux densities for SN1957D were obtained at each frequency using an elliptical aperture 
centered on SN1957D, with major and minor axes equal to twice the restoring beam FWHM 
and rotated to the same position angle as the restoring beam. The flux densities are 230±10 
/iJy and 190±10 /iJy at 5.5 GHz and 9 GHz respectively, giving a spectral index a = 
— 0.4±0.1 (F v oc u a ). This spectral index is more typical of SNRs in which the radio emission 
i s dominated by shocks than those in which the radio emission is dominated by a pulsar 
( lGreenll2009i ). There are no previous radio measurements of the flux from SN1957D at 9 
GHz, and therefore we cannot directly compare radio spectral indices from these observations 
to those made by others at 6 and 20 cm, the most recent measurement of which i s a = 
-0.26 ± 0.14 based on data taken in 1998 iMaddox et al.lbood : Istockdale et~alll2006[ ). 



The time evolution of the radio flux from SN1957D is shown in Fig. [3J A broken power 
law of the form 



F v OC t t < ^break 

52 t>ty 



F v oc t 



" break 



reproduces the data, and captures the fact that the rate of declined has steepened over time. 
The power-law slope after the break is well constrained (<5 2 = 3.97+q;2q), but the location 
of the break (tbreak < 33 years since the SN explosion), and the power-law slope before the 
break (8% < 1.5) are not as well constrained. In Fig. [3] we have assumed tbreak = 33 years 
since the SN explosion — the upper limit for tbreak- However, we might expect th at the break 
is coi ncident with the rapid drop in [O III] flux seen in the period 1988-1989 ( ILong et al. 
19921 . see also Section 0] and Fig. [3j bottom panel), in which case the initial decline would 
be much flatter than the model fit in Fig. [3j with slope of approximately 5\ ~ 0.7. 



4. Optical Observations — Reductions and Analysis 
4.1. Imaging Observations 



Long et al.l (119921 ) discussed the early time-evolution of the [O III] emission from SN1957D, 
noting that the flux had declined by a factor of five between 1987 and 1991. We have contin- 
ued to obtain narrow-band [O III] images of the SN over the years to monitor its evolution. 
An observing log for the various datasets we have available is presented in Table [TJ We 
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have reduced all of the ground-based data using standard IRAF0 procedures, including flux 
calibration based on observations of standard stars from the list of iHamuy et al. ( 1992 ) and 



the known transmission functions for the filters. We subtracted scaled continuum images 
from the emission-line ones and carried out conventional aperture photometry to measure 
the fluxes. In the upper right panel of Fig. [TJ we show a small portion of one of these data 
sets, from the 6.5m Magellan/IMACS continuum-subtracted emission-line images. 

The region of M83 containing SN1957D was observed with the HST WFC3 in 2010 
March as part of the Early Release Science program after Servicing Mission 4 see Table 
[2]). We obtained both narrow-band emission-line and broadband continuum data from the 
Mikulski Archive for Space Telescopes (MAST) at STScI and, starting with data redu ced 

to 
A 



through the standard WFC3 pipeline, used MultiDrizzle ( IFruchter h Sosey et al.l 12009 
combine the various exposures and put all the images on a common coordinate system 
scaled version of the F547M (narrow V) continuum image was used to remove most of the 
stellar emission from the narrow-band images to produce pure emission-line images in Ha, 
[S II], and [O III]. These data were used to produce the lower right panel in Fig. [TJ which 
is directly comparable to the Magellan data presented above. The HST data clearly resolve 
the point-like [O III] SNR emission from the surrounding more diffuse gas, and in particular 
a faint arc of Ha emission to the south and east of the SNR position. We have included 
our photometric measurements for the [O III] 5007 A emission along with the ground-based 
measurements in Table [TJ 

The HST broadband images of the region are also very interesting. In the lower left 
panel of Fig. [TJ we show an aligned 3-color representation of the HST U-B-narrow V data, 
which demonstrates the presence of a very blue association of stars just to the northeast 
of the SNR position. Furthermore, careful alignment of the continuum and emission-line 
images shows that exactly at the position of the SNR is a remaining star, whose color is 
somewhat redder than most of the other stars in the adjacent association. This figure also 
shows that the extinction in the region is patchy and variable, extending right to the edge 
of the association on the west side, and possibly overlying the SNR position itself. We do 
not know if this star is associated with the SN precursor (e.g. a possible companion star) or 
whether it is a chance alignment along the line of sight, but its position is coincident with 
the SNR position within the measurement error. The association will be discussed further 
below. 



5 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by AURA, Inc., 
under cooperative agreement with the National Science Foundation 



6 These are the same data used by iMilisavlievic et al.l (|2012f) . but we have done the reductions and pho- 
tometry independently. 
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The precision of the flux values in Table [T] is limited primarily by our ability to dis- 
tinguish SN1957D from the nearby H II-region emission, and to accurately measure the 
underlying background. Both problems are, of course, worse under conditions of poor see- 
ing, but even with the HST WFC3 data there is still some background contamination. In 
order to minimize the confusion problem, we have used minimum aperture sizes (diameter 
FWHM) for the SN1957D photometry and have made aperture corrections derived from 
data on several isolated bright stars in each image. We have also experimented with different 
choices for background regions, and include variation from these, in addition to purely statis- 
tical errors, in the uncertainty estimates for Table [TJ As noted in the table, there are minor 
differences in the pass bands for the different [O III] filters used for different observations, 
but all pass most of the broad 5007 A emission, and exclude most of the 4959 A line. All 
the flux values also include the narrow component of the 5007 A line, since there is no way 
to distinguish between broad and narrow components in our images. 

The time evolution of the [O III] flux derived from the above observations is shown in 
the bottom panel of Fig. [3j All of the measurements are consistent with an [O III] flux of 
about 10~ 15 ergs cm -2 s _1 between 2002 and 2011, about 2.5 times smaller than the flux 
measured in 1991. Clearly the rapid decline in the flux seen in the earlier data has halted. 
Unfortunately we do not have any images between 1991 and 2002, but interpolating the 
existing data, it appears consistent that the rapid decline had essentially leveled-out by the 
early 1990's. 



4.2. Spectroscopic Observations 

As part of a Gemini observing program 2011A-0436, Winkler PI) to obtain spec- 
tra of SNRs in M83, we obtained a spectrum of SN1957D using the Gemini Multi-Object 
Spectrograph (GMOS) on the 8.2m Gemini South telescope on 2011 April 8 (UT). Our ob- 
servations consisted of three 40-minute exposures using the 600 line blue grating (G5323) in 
multi-slit mode. For SN1957D, the custom slit was 1"25 wide and 10" long, oriented E-W. 
The spectrum, extracted from a 1.5" region along the slit and covering the wavelength range 
4700 - 7300 A, is shown in Fig. HI The original two-dimensional spectrum has a blue con- 
tinuum coincident with SN1957D but extending extending well beyond it to the east, which 
we attribute to the blue stellar association shown in Fig. [T] (lower left panel); for clarity we 
have fitted and subtracted that from the spectrum in Fig. HI The spectrum shows broad 
lines from [O III] AA4959,5007, [O I] AA6300,6363, [S II] AA6717,6731, and [Ar III] A7135, 
as well as narrow lines from Ha, H/?, [He I] A5875, [O III] AA4959, 5007, [N II] AA6548, 6583, 
[S II] AA6717, 6731, and [Ar III] A7135, as listed in Table 3. Simultaneous fits for the broad 
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components of the [O III], [O I], and [S II] lines gives a velocity width of 2730 ± 100 kms -1 
(FWHM). From the HST and ground-based imagery available, we conclude that some of the 
narrow-line emission probably arises from the near-coincident diffuse emission located < 1" 
to the east and south of SN1957D, since the two cannot be fully resolved at ground-based 
resolution. However, it appears that some of the narrow-line emission does indeed stem from 
SN1957D itself. We attribute the fact that the integrated [O III] flux values (Table 3) are 
lower than the those from imaging (Table 1) to the 1725 slit's having excluded a portion 
of the SN1957D flux, and/or to a slight mis-positioning of the slit. The flux values from 
spectroscopy shown in Figj3] are systematically lower than those from imaging at similar 
epochs, probably for the same reason. 

The broad-line profile appears to be non-Gaussian, probably reflecting the distribution 
of shocked ejecta within the young remnant. Fig. shows the [O III] lines at various epochs 
from 1989 through 2011. In addition to the 2011 GMOS s pectrum, these include spectra 
from epochs 1989 and 1991 reported by lLong et al.l (119921 ). and a previously unreported 
spectrum obtained from the 4m Blanco telescope at CTIO in 1998 June. The broad line 
profile appears qualitatively rathe r similar at all these epochs, and indeed at all epochs 



since SN1957D's optical recov ery (lLong et al.l Il989l ; iTuratto et al.l Il989l ; ICappellaro et al. 



19951 ; iMilisavljevic et al.l 120121 ). The combined flux in the broad component of the [O III] 
AA4959, 5007 lines in 2011 is 4.2 x 10~ 16 ergs cm -2 s -1 ; all of the broad and narrow line fluxes 
are given in Table El 



4.3. The Stellar Population near SN1957D 

The stars in the region of M83 that contains SN1957D are very blue, indicative of a 
young association (see Fig. [I] and the expanded view in Fig. [6]). There is a fair amount of 
dust in the area, but the very blue colors of most of the stars suggest that the association 
is mostly on a fairly clear line of sight with minimal extinction. There are some red giants 
in the general vicinity (the bright orange stars in Fig. [HD, but none is within 2'.'2 (50 pc) 
of the SNR and none is within the concentration of stars that defines the association. The 
absence of red giants suggests a relatively young age for the association. To quantify this, 
we carried out DAOPHOT aperture photometry of stars in the association, as defined by a 
rectangular box with dimensions 52 x 35 pixels (2'.'1 x l'/4 or 46 x 31 pc), with the SNR near 
the right side of the rectangle (i.e., the upper left two-thirds of the field-of-view in Fig. 
We used an aperture radius of 3 pixels (0'.'12), background sky annulus from 10 to 13 pixels, 
and the standard encircled-energy expression to account for the missing flux outside the 3 
pixel aperture. 
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We corrected the magnitudes for Galactic foreground extinction using a value of Ay = 
0.210 and produced a color-color diagram shown in the left panel in Fig. For comparison, 
we also plot the Padova stellar isochrones for solar metallicity stars of age 4 and 10 Myr 



( iBertelh et al.l 11994 , 120091 ) . The stars in the M83 stellar association are near the blue tip of 
the stellar isochrones, nearly coincident in the color-color diagram, as appropriate for very 
young stars. There is more scatter in the V — I direction than the U — B direction since 
these are very blue stars (i.e., some are quite faint in the I band). Only stars brighter than 
My = -4 are included in our analysis. The median value of V — I is -0.2, which matches the 
models very well, suggests that there is minimal internal extinction for most of the stars. 
However, the six reddest stars (solid square symbols in Fig. [7J do appear to be aligned 
roughly along the vector expected from reddening (dashed line). Hence we have correcte d 



these stars individually for extinction using the technique described in iKim et al.l ( 120121 ). 
The corrected data are shown by the open squares. 

The right panel of Fig. [7J shows the corresponding color-magnitude diagram for the stars 
in the SN1957D association, using the same symbol definitions as in the left panel. The 
locations in the diagrams are shown both uncorrected and corrected for internal reddening. 
After the data for the six reddest stars have been corrected for reddening, essentially all of 
the data are consistent with the stars being on the main sequence of a 4 Myr isochrone. A 
10 Myr isochrone is also included, representing the oldest possible age for the association 
(i.e., the brightest star would have transitioned to a red giant). The turnoff masses for the 
4 and 10 Myr isochrons are 38 and 17 M Q respectively, both consistent with that expected 
for a core-collapse SN. 



Discussion 



SN1957D is not the only decades-old SN (or very young SNR) to have detectable X-ray 
emission. As a result of the increasing sensitivity of X-ray telescopes, primarily Chandra 
and XMM-Newton, about a dozen old SNe have been detected as X-ray sources a decade or 
more after the SN explosion ( jSoria fc Pernall2008l ; iDwarkadas fc Gruszkol 120121 ) . The oldest 
extragalactic SN that has been detec ted in X-rays today is SN1941C at a luminosity of 
5 x 10 37 ergs s _1 (jSoria fc Pernal 120081 ) . At age 24 years, SN1987A has an X-ray luminosity 



of about 4 x 10 ergs s , though unlike SN1957D, its spectrum is dominated by soft X-rays 



in the range 0.3-2 keV that arise from a shock-heated thermal plasma (see, e. g. iPark et al. 



7 T he average of estimates of the Galactic extinction by 
dl998h . 



Burstein fc Heiles 



(1982) and 



Schlegel et al, 
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2 111 , and references therein). Others range in X-ray luminosity from about 10 37 ergs s 1 



to 10 41 ergs s 1 . Most of these have X-ray spectra th at are soft, and in cases where the 



sourc es are sufficiently bright, such as SN1979C in M 100 (jlmmler et al.ll2005t iPatnaude et al. 



201 ll ) and SN1993J in M81 ( jChandra et al.l 120091 ) , evidence of emission from lines. The 
fluxes are usually constant or declining with time (IDwarkadas fc Gruszkoll2012l ). Their X- 
ray properties are normally interpreted as the interaction of the SN shock with circumstellar 
material, although IPatnaude et al.l ( 1201 ll ) have suggested that SN1979C has a non-thermal 
component as well. All are assumed to be core-collapse SNe. The lower limit of about 
10 ergs s _1 for a detected X-ray SNR reflects, for the most part, the sensitivity limits of 
the existing observations given the distances of nearby galaxies. X- ray emission has no t been 
reliably detected from any Type la SN in this age range (see, e.g. iHughes et al.l 120071 ). 



The X-ray spectrum of SN1957D is featureless and hard, unlike the other SNRs in this 
age range that have been observed in X-rays and unlike the majority of Galactic SNRs. 
Objects that do resemble SN1957D are pulsar-powered SNRs such as the Crab Nebula and 
the SNR 0540-69.3 in the Large Magellanic Cloud. The Crab Nebula is 958 years old; 
SNR 0540-69.3 i s tho ught to be between 700 years old from the expansion o f the filaments 



dWil 



iams et al.l 120081 ) and 1700 years old from the pulsar spin-down age ( jGradari et al. 



20111). X-rays from bo th are do minated by non-therm al emission with photon indices of 2.1 



( IWeisskopf et al.ll2010l ) and 2.0 ( IHirayama et al.ll2002f ). respectively, compared to SN1957D 
with 1.4+Qg . The 0.3-8 keV luminosities of all three objects are about the same, ~ 
1 — 2 x 10 37 ergss -1 . Thus, it is plausible that X-ray emission in SN1957D is powered 
by a pulsar. 



The column density that we derive from the X-ray spectrum of SN1957D is quite high, 



2.01];' ^ x 10 22 cm 1 1 m uch hi gher than expected from for eground absorption in the Galaxy 



or in M83. As discussed by iFransson fc Chevalierl (119871 ) for SN1987A, cold ejecta from a 
core-collapse SN explosion provide substantial opacity for many yea rs after the explosion. 
Indeed , if the event that generated SN1957D is similar to that used by lFransson &: Chevalier 
( 119871 ) for modeling SN1987A, one would predict that today the energy with optical depth 
unity would be roughly 5-6 keV. Thus, if X-ray emission in SN1957D is powered by a pulsar 
it is not surprising that the source would be highly absorbed. 

The points marked with an X in the color-color and color-magnitude diagrams shown 
Fig. [7] represent the star that is coincident with the SNR position. With a reddening A^ of 
about 2.1, this star is the reddest of the stars in the association. The other five stars we 
dereddened to place on the stellar isochrones have an average Ay of 1.4, with a maximum 
value of about 1.8. As mentioned earlier, SN1957D is on the edge of the association and 
local extinction variations could account for the higher reddening of the star. A more exotic 
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idea would be that this is the companion to the SN and that localized dust formed in the SN 
ejecta is affecting the reddening. A reddening of Ay ~2.1 co rresponds to a co l umn d ensity 
of 3.3 x 10 21 cm~ 2 if the gas properties are typical of the ISM lSavage &: Mathisl ( 119791 ). This 
is far less than the column density, 2.0!^ x 10 22 cm -2 , we derive for the X-ray source, but 
there is no reason to believe that the properties of the gas in SN ejecta resemble that of the 



ISM. Indeed the dust content of young SNe have turn ed out to be c 



predictions based on normal gas-to-dust ratios (see, e.g. iNozawa et al.ll2010l ). As a result, we 



uite low compared to 



have no observational way to distinguish between the possibility that the star is physically 
associated with the SN, or that its higher reddening is due to local extinction variations. 

Anomalous X-ray puls ars ( AXPs) , thou ght to be ultramagnetized (B> 3 x 10 14 G) iso- 



lated neutron stars (see, e. g. lMereghettill2008l . for a review of their properties), are in principle 
another way to power emission from young SNRs . They are known to be associated with 
young SNRs, and, since normal pulsars, with lower magnetic fields, generate pulsar wind 
nebulae, one suspects that they should be able to produce pulsar wind nebulae as well. 
However, the average luminosity of AXPs is typically 10 34-36 ergs s _1 , less than is needed to 
power the observed emission from the remnant of SN 1957D and AXP X-ray spectra usually 
show steeper spectral indices than we observe from the remnant of SN1957D. Furthermore, 
at present, the exis tence of (even faint) pulsar wind nebulae around AXPs is controversial 
ffYbunes et alJl2012h . While we cannot rule out an AXP explanation for what we observe in 
SN1957D, given the state of our knowledge of AXPs, there is no need to invoke this expla- 
nation. The only obvious way to distinguish between a normal pulsar and an AXP would 
be to detect it directly, most likely by observing a soft 7 ray outburst in the direction of 
SN1957D. This is not possible with today's instrumentation. In the absence evidence for an 
AXP in SN1957D, a more normal pulsar/pulsar wind nebula, similar to the one that exists 
in the Crab Nebula is favored. 



Patnaude et al.l (120111 ). in order to explain a non-thermal component in the X-ray spec- 
trum of SN1979C, discuss the possibility that the emission might arise from a black hole 
accreting material from a SN fall-back disk or from a binary companion. The X-ray lu- 
minosity of SN1979C is 6.5 x 10 38 ergss -1 , comparable to the Eddington luminosity for a 
stellar BH, which was one of their reasons for arguing that this system might require a black 
hole. SN1957D is 30 times less bright in X-rays and so a rapidly spinning pulsar of the type 
known to exist in the Crab Nebula and SNR 0540-69.3 can provide the needed luminosity for 
thousands of years. A black hole origin for the X-ray emission in SN957D cannot be ruled 
out of course, but again, there is no indication that it is required in SN1957D, and therefore 
a more normal pulsar wind/nebula is favored for the X-ray emission. 



In contrast to its X-ray properties, the radio and optical properties of SN1957D are not 
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unusual for very young core-collapse SNRs. Radio emission in young core-collapse SNRs is 
generally interpreted as synchro tron emission ge nerated as the primary SN shock interacts 
with the circumstellar medium (j Chevalier 1982 ). The radio luminosity of SN1957D is not 
exceptional. A secular decline in radio flux for a core-collapse SN as F v oc t -0,7 is expected if 
the circumstell ar medium is th e stellar wind from a supergiant progenitor, where the density 
scales as r~ 2 (j Chevalier! 1 19821 ; IWeiler et al.ll2002l ). The decline should steepen when the 
blast wave overruns the edge of the circumstellar material, and then flatten as it enters the 
undisturbed ISM. Break s in the radio light curve have been meas ured for the recent type lib 
SNe 1993J and 2001gd JWeiler et alil2007l : Istockdale et alil2007h . In both cases, the breaks 
are attributed to a decrease in the density of the circumstellar medium associated with a 
sudden increase in the mass-loss rate in the years prior to the explosion. The abrupt change 
in the radio flux is particularly well documented for SN1993J, beginning 8.5 yrs after the 
explosion. This appears to be what has happened in the case of SN1 957D. Between 1984 
and 1990, it declined roughly at the canonical rate ( jCowan et al.lll994l ). However based on 
our new ATCA data, the rate has been declining as t~ 4 since about 1990. T his suggests, 



assum ing a wind speed of 10km s 1 and a velocity for the shock of 15,000km s 1 ( IWeiler et al. 



20071 ). that the increase in mass loss from the progenitor of SN1957D occurred about 8000 
yrs prior to explosion. 

Similarly, as noted in Section (TJ the optical spectra of young core colla pse SNRs are often 



but n ot always dominated by emission from forbidden lines of oxygen ( IMilisavljevic et al. 



20121 ). The optical emission from young core-collapse SNRs, whether dominated by oxygen 
lines or not, is understood to arise from shock interactions within the eject a and between the 
forward shock and the circumstellar medium (j Chevalier fc Franssonlll994l ). For the oxygen- 
dominated remnants, such as SN1957D, the best picture to date one that involves the inter- 
action between fast-moving, dense oxygen-dominated fragments from the explosion and less 
dense, perhaps more evenly distributed SN ejecta of the reverse shock ISutherland fc Dopita 
(119951 ). (This reproduces the observed line ratios better than simpler models involving the 
reverse shock propagating back into the ejecta.) The interaction between the fragments and 
the more tenuous material creates a bow-shock around the fragments, or knots, and drives 
shocks into the fragments themselves. Optical emission arises from the shocked oxygen- 
rich ejecta in the knots. Optical spectra of SN1957D are similar to those of other O-rich 
SNRs, including Cas A, G292.0+1.8 and E0102-72.3, that lack luminous pulsar/pulsar wind 
nebulae. 

Nevertheless, in view of the unusual X-ray spectrum of SN1957D and the suggestion 
that its X-ray emission arises from a pulsar/pulsar wind nebula, it is worthwhile asking 
whether the radio or optical emission could be due to a pulsar as well. 
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The radio flux from the pulsar wind nebula within SNR 0540-69.3 is 520 mJy at 5 
GHz ([Manchester et al.lll993l ). Scaled to the distance of M83, this flux would be 61 fiJj, 
co mpared to the 230 //Jy at 5 GHz that we measure from SN1957D. Using parameters given 
by iGreenl (120091 ). the Crab Nebula would have a flux of about 120 /j,Jj at 5 GHz at the 
distance of M83, about half that of SN1957D in the latest ATCA observations. Thus, while 
the radio luminosity of SN1957D is larger than either of these two objects, it is not so large 
as to make implausible the possibility that the radio emission is powered by a pulsar (if 
we include the pulsar wind nebula in the estimate of the luminosity). However, a problem 
with the hypothesis that a pulsar dominates at radio wavelengths is how to reconcile in a 
natural manner the time evolution of the flux at radio and X-ray wavelengths. We know, 
based on the upper limits to X-ray flux from the 2000-2001 Chandra observations, the X-ray 
luminosity in 2000-2001 was about the same as or possibly lower than in 2011. Although 
synchrotron lifetimes for electrons differ at different wavelengths, very special circumstances 
would appear to be required to create a rapidly declining flux at radio wavelengths and a 
near constant (or increasing) flux at X-ray wavelengths. 

The optical lines generated in pulsar-dominated supernova remnants, such as the Crab 
Nebula, a re produced primarily b y pho toionization of ejecta (and the ISM) by synchrotron 



emission. iChevalier fc Franssonl (119921 ) discussed the development of UV-optical-IR emis- 
sion lines in pulsar-dominated core-collapse SNR. In very young SNRs, the lines are nar- 
row and permitted lines are favored as the ionization progresses into the expanding ejecta, 
but the velocity widths broaden with time and forbidden lines, especially [O III] become 
more prominent as the density drops. The optical spectra of the Crab Nebula and SNR 
0540-69.3 are similar to that observed in SN1957D, with broad emission from the lines of 



oxygen and sulfur (see, Fig. 6 of iMorse et al.ll2006l . for a comparison), although Ha+[N II] 
AA6548, 6583 is much more obvious in the integrated optical spectrum of the Crab. The 
[O III] AA4959, 5007 emission-line fluxes in SN1957D have been nearly constant since 2002, 
which is consistent with the X-ray fluxes, but if we speculate that the [O III] AA4959, 5007 
flux scales with emission from the pulsar, then the pulsar/pulsar wind nebula must have 
bee n dramatically brighter prio r to 1990, which is not easily understandable. Moreover, in 
the lChevalier fc Franssonl (119921 ) models, [O III] AA4959, 5007 should brighten with time and 
the line widths should broaden with time. That is not what is observed in SN1957D. 

A final possibility is that the radio emission, optical emission, and the X-ray emission 
all arise from the shocks between ejecta and circumstellar material. But in this case, one 
would again need to explain the hard, featureless spectrum and the fact that the X-ray object 
was not brighter in the 2000-2001 Chandra observations. Unfortunately, it is unlikely that 
a better X-ray spectrum will be acquired any time in the near future, though it should be 
possible to monitor the X-ray flux. Possibly the X-ray (and radio emission) arise primarily 
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from synchrotron emission at the shock front. However, to our knowledge, none of the SNRs 
of this type is nearly as bright as SN1957D in X-rays and none has emission lines arising 
from ejecta. Perhaps the most interesting of these objects is the Milky Way SNR Gl.9+0.3, 
which, with an estimated age of 140±30 years, is the youngest known SNR in the Galaxy. It 
has a very hard non-thermal spectrum, similar to that observed in SN1957D, but its X-ray 
luminosity is ~ 3 x 10 34 ergs s _1 , about 50 times fainter than S N1957D, consistent with its 



being the remnant of a Type la explosion ( {Reynolds et al.ll2009i ). The alternative, emission 



from a very hot thermal plasma, also seems hard to justify, given the X-ray spectra of other 
young core-collapse SNe. 

In conclusion, while we caution the reader to keep in mind the fact that the X-ray 
spectrum contains only 140 counts, the existing data favor a picture of SN1957D in which 
the X-ray emission arises from a pulsar/pulsar- wind nebula, while the optical and radio 
emission arise primarily in shocks in the ejecta and circumstellar gas. In this picture, X- 
rays are produced in a region within the cloud of cold ejecta and are thus heavily absorbed. 
The optical emission arises in the vicinity of the reverse shock, but outside the cold, freely 
expanding ejecta. The radio emission is produced at the shock front. 



6. Summary 

We have detected X-ray emission from the remnant of SN1957D for the first time using 
a very deep Chandra observation of M83. The X-ray spectrum of the source is unusually 
hard and has unusually high column density for a young SNR. The X-ray spectra suggest 
that the SNR is dominated by an X-ray pulsar with substantial absorption in the ejecta. 
The SN was almost certainly the core-collapse explosion of a massive progenitor, based on 
its association with a young star cluster with an age of 4-10 million years, the high star 
formation rate in the spiral arms of M83 generally. At this age, the turnoff mass is at 
least 17 M Q , exactly where one would expect a core-collapse SN. Optically, the emission is 
dominated by [O III] AA4959,5007, suggesting a shock interacting with ejecta that is highly 
enriched in oxygen. The shape of the optical spectrum has not evolved significantly since the 
first spectrum was obtained in 1987, 30 years after the explosion, but the flux has declined 
at both optical and radio wavelengths indicating that the density of the medium into which 
the shock is expanding has also declined. Given the unusual nature of the X-ray spectrum, it 
is important to continue to monitor the evolution of this very young SNR at all wavelengths. 
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Table 1. SN1957D Imaging Observations and Fluxes 



Date (U.T.) 


Telescope 


CCD 


Line 


Filter 

T / S \ 

A c (A) 


AA a (A) 


Exposure (s) 


Observers 


[OIII] 5007 A Flux 
10 -16 ergs cm~ 2 s~ 1 


1987 Apr 26 a 


LCO 2.5 m 


Tl 


[0 III] 


5020 


54 


2 X 1000 


KSL, WPB 


120 








Continuum 




1 nn 

1UU 




W. Krzcminski 




1991 Apr 18 b 


CTIO 4.0 m 


TEK IK #1 


[O III] 


5020 


51 


3 x 600 


PFW, KSL 


21 








Continuum 


4770 


100 


3x500 






2002 Mar 20-23 


CTIO 0.9 m 


TEK 2K #3 


[O III] 


5006 


60 


4 x 1000 


PFW, KSL, 


8.9 ±2.6 








Continuum 


5135 


90 


4 x 500 


C. Reith 




2007 Apr 15 


SOAR 4.1 m 


SOI 


[O III] 


5006 


60 


4 x 600 


PFW, KSL 


8.1 ± 1.6 








Continuum 


5316 


161 


7 x 200 






2009 Apr 26 


Magellan I 6.5 m 


IMACS 


[O III] 


5007 


50 


7 x 600 


PFW, KSL 


8.0 ± 1.0 








Continuum 


5316 


161 


7 x 200 






2010 Mar 19 


HST 


WFC3 


[O III] 


5010 


65 


3 x 828 


R. O'Connell 


8.9 ±0.5 



HLong et al.l lll989h 
^Long et al.1 il99j ) 
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Table 2. HST Image Summary" 



Band 


Association ID 


Filter 


Exposures 


Total Exposure 

00 


U 


IB6WB1060 


F336W 


3 


2560 


B 


IB6WB1050 


F438W 


3 


1800 


V 


IB6WB2070 


F547M 


4 


1203 


I 


IB6WB2060 


F814W 


3 


1213 


[OIII] 


IB6WB2050 


F502N 


3 


2484 


Ha 


IB6WB2010 


F657N 


4 


1484 


[SII] 


IB6WB3030 


F673N 


3 


1770 



a Imagcs obtained in 2010 March as part of Early Release Observations 
for WFC3. Program ID: 11360, PI: Robert O'Connell 
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Table 3. SN1957D Emission-line Fluxes, 2011 April 



Flux b 

Ion A (A) a 10~ 17 ergs cm -2 s _1 



Broad Lines, AA = 2730 kms" 1 



[0 III] 


4959 


12.1 


[0 III] 


5007 


30.2 


[0 I] 


6300 


10.1 


[0 I] 


6363 


3.0 


[S II] 


6716 


8.5 


[S II] 


6731 


5.9 


[Ar III] 


7135 


(5.7) 


Narrow Lines (unresolved) 


H/3 


4861 


14.8 


[0 III] 


4959 


4.7 


[0 III] 


5007 


11.8 


[He I] 


5875 


3.2 


[0 I] 


6300 


(1.0) 


[N II] 


6548 


12.4 


Ha 


6563 


76.2 


[N II] 


6583 


37.4 


[S II] 


6716 


9.2 


[S II] 


6731 


6.9 


[Ar III] 


7135 


3.3 



a Rest wavelength 

b Values in parentheses are highly uncer- 
tain. 
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Fig. 1. — A four panel context figure showing the region in the immediate vicinity of 
SN1957D. At upper left is a three-color Chandra image, where red is 0.35 - 1 keV, green is 1 
- 2 keV, and blue is 2 - 8 keV. At upper right is a three-color Magellan optical continuum- 
subtracted emission line image, taken in exceptional Of! 4 seeing, where red is Ha, green is 
[S II], and blue is [O III]. At lower right is the same as the upper right panel, but for the 
HST WFC3 emission line data. The lower left panel is the HST WFC3 continuum plus Ha 
image of the same region, provided by Zolt Levay (STScI), where red is F547M, green is B 
band, and blue is U-band, with Ha also shown as red/pink. (Note the star subtraction on 
the [S II] data in the lower right panel is imperfect, leaving some stellar residuals that can be 
used for reference in comparison with the lower left panel.) The green circles indicate normal 
ISM-dominated SNRs in the vicinity, and the yellow circle indicates the SN1957D position. 
The hard X-ray source is positionally coincident with the optical SNR which is dominated 
by [O III]; the star indicated by the arrow in the lower left panel is also coincident with the 
optical SNR (see also Fig. |6]for an expanded view of this region). North is up and east is to 
the left, and a scale bar is shown in the upper left panel. 
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Fig. 2. — X-ray spectrum of SN1957D fitted with an absorbed power law. Datapoints and 
C-statistics residuals are shown. After allowing for Galactic absorption of 4 x 10 
find a photon power law index 7 of 1.4+^ and A^(M83) of 2.0±{% 
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Fig. 3. — The time evolution of the radio flux (upper panel) and the optical [O III] A5007 
emission-line flux (lower panel) of SN1957D. The horizontal axes are identical and the verti- 
cal axes have the same dynamic range in the two panels. In the optical panel, filled symbols 
represent flux values from imaging; open symbols represent values from spectra, which typ- 
ically miss some flux. 



-26 - 




5000 5500 6000 6500 7000 

Wavelength (A) 

Fig. 4.— Optical spectrum of SN1957D, as observed in 2011 April with GMOS on Gemini 
South. The remnant shows broad lines from [O III] AA4959, 5007, [O I] AA6300, 6363, and 
[S II] AA6717, 6731. Narrow emission lines are also seen, but they most likely arise from 
photoionized gas seen to the east of the SNR in the HST images. 
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Fig. 5. — The emission- line profile in [O III] AA4959, 5007 region of the spectrum a s observed 
in 2011 April with GMOS (Fig. gj, and from 1989, 1991 fboth lLong et al.lll992h and from 
1998 (previously unpublished CTIO Blanco data). The broad line emission has decreased 
over time (note scaling), but the profile has remained qualitatively similar. The narrow lines 
were largely removed from the 1989 and 1991 spectra, but have been retained (but truncated) 
in those from 1998 and 2011, where the higher resolution minimizes broad-narrow confusion. 




Fig. 6. — The stellar association in which SN1957D is located as observed with HST. The left 
panel is a 3-color composite of the emission line gas, where red is Ha, green is [S II], and blue 
is [O III]. The right panel is a portion of the continuum plus Ha image of the same region 
described more fully in the caption of Fig. [U Stars within the rectangular shaped region 
were used to establish the age of the stellar association along the line of sight to SN1957D. 
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Fig. 7. — Color-color (left panel) and color-magnitude (right panel) diagrams for stars in 
the stellar association along the line of sight to SN1957D. The positions of the stars in the 
assocation are compared t o stellar isochrones for solar metallicity stars of age 4 and 10 Myr 
(IBertelli et al.lll994l 120091 ). For stars with significant reddening, both the uncorrected and 
reddening-corrected positions are shown. The stellar association is likely to be as young as 
4 Myr and is almost certainly less than 10 Myr. The X indicates the object coincident with 
SN1957D. 



